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Abstract 

Summations and relations involving the Hurwitz and Riemann ^- 
functions are extended first to Barnes (^-functions and then to 
functions of general type. The analysis is motivated by the eval- 
uation of determinants on spheres which are treated both by a di- 
rect expansion method and by regularised sums. Comments on ex- 
isting calculations are made. It is suggested that the combination 
C'{—n) + H^({—n), where is a harmonic number, should be taken 
as more relevant than just n). This leads to a Kaluza-Klein 
technique, providing a determinant interpretation of the Glaisher- 
Kinkelin-Bendersky constants which are then generalised to arbitrary 
^-functions. This technique allows an improved treatment of sphere 
determinants. 
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1. Introduction 



The Barnes ^-function, [1], has enjoyed only sporadic mathematical interest 
since it was first introduced in 1903, [2-4]. It generalises the Hurwitz C-fmiction, 
has number theory applications and involves a natural generalisation of the Euler 
F-function, examples of which had already appeared, [5-7] and were developed more 
directly by Bendersky [8]. 

Some of these works discuss regularised products, [9-12], and the related Lapla- 
cian determinants, [13-15], on various spaces from varying points of view. Explicit 
results are usually obtained only when the eigenvalue spectrum is known or suffi- 
ciently powerful geometrical or analytical information is available. 

As might be expected, particularly detailed investigations have been made 
in the case of spheres. Explicit d-sphere determinant expressions are given, and 
plotted, in [16], obtained by a direct method, ancillary to a discussion of the orbifold 
quotient, S^/T. Choi [17] analyses the 3-sphere using a factorisation technique 
founded on the method in Voros' basic paper [14] while Choi and Srivastava [18] 
improve this calculation. Quine and Choi [19], with a systematic method, give 
explicit forms for the full d-sphere. Kumagai [20] attempts to give a corrected 
version of Vardi's analysis, [13], for the (i-sphere but does not calculate beyond the 
4-sphere."^ Further remarks are made later. 

During the course of these investigations, identities and relations appear which 
wc would like to systematise in a certain way and enlarge upon. Some considerable 
amount of work in this area has appeared in the mathematical literature in the past 
few years and we wish to draw attention to relevant work by physicists which is 
often overlooked and which might be useful. 

Applications of the Barnes function in physics embrace the Casimir efi^ect 
around cosmic strings [21], on orbifolded spheres, [22,23], integrable field theo- 
ries [24-27], and statistical mechanics, [28]. The Barnes ^-function also arises for 
the higher-dimensional harmonic oscillator and is useful in connection with Bose- 
Einstein condensation and trapping [29,30]. 

As noted in [16], it is clear there is some error in Vardi's manipulations. This can be traced 
to two simple oversights on p. 504 of [13]. The first is in the proof that there are no logfc 
terms. The upper limit on the k sum cannot be changed from n — 1 to n — 2 because there is a 
contribution when d = 0. The second slip is an incorrect interchange of the d and r summations. 
Vardi produces a general formula for the d-sphere determinant in terms of the derivatives of 
the Riemann (^-function at negative integers, cf [16,19], which he then converts into multiple 
r-functions, r„(l/2). 
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The structure of this paper is as follows. The Barnes ^-function is introduced 
via scalar ^-functions on spheres. The determinant is next looked at, which focuses 
attention on certain constructs and relations. These further motivate an investiga- 
tion of the 5'eneroi ^-function, ^ A~^, and we then link up with the regularised prod- 
uct approach giving a Kaluza- Klein interpretation of generalised Glaisher-Kinkelin- 
Bendersky constants. The paper proceeds as a series of generalisations. 

Determinants appear frequently in field and string theory. No attempt will be 
made to justify their computation nor to detail their general history. For the 2- 
sphere the first computations were by Hortacsu et al [31] and by Weisberger [32,33]. 
The topic of C-functions on spheres, and symmetric spaces in general, also has a 
long record. (See Camporesi, [34], for a useful survey and results.) 



2. The sphere C— function and the Barnes function. 

The general definition of the Barnes (^-function is, 

r (. n\^) r exp(-ar)(-T)--^ 

-^L Ui^i {1 - eM-diT)) 

oo ^ 

= > 7 TT- , Re s > d, 

m=0 ^ ' 

where we refer to the components, dj, of the d- vector, d, as the degrees or quasi- 
periods. For simplicity, we assume that the d^ are real and positive. If a is zero, the 
origin, m = 0, is to be excluded. The contour, L, is the standard Hankel one. 
On the (i-sphere, consider the Laplacian~type operator, 

-A + CR-a'', (2) 

where ^ = {d—l)/4:d and a is a parameter introduced for convenience of discussing 
several specific cases at once. The value a = corresponds to conformal coupling in 
d+1 dimensions, a = 1/2 to conformal coupling in d-dimensions and a = (d — l)/2 
to minimal coupling i.e. just the operator —A. In the first case the eigenvalues are 
perfect squares, 

^{m + d-2)^, m=l,3, .... (3) 

It is best to think of the mode set on the full sphere as the union of Dirichlet and 
Neumann mode sets on a hemisphere, despite the apparent extra complication. The 



2 



reason is that these individual ^-functions are Barnes functions, with all degrees 
equal to unity, d = 1, as shown in [22]. Specifically, for a = 0, 

C^(.) = Crf(2., (d-l)/2|d), 

C^(.) = C,(2s,Erf.-(rf-i)/2|d), 

where the parameters in the ^-functions have been left general because, although 
for the full-sphere, and hemisphere, we need only unit degrees, we wish to retain 
the general case for as long as possible as it applies for the other orbifold factorings 
of the sphere^. 

If the sum over the vector m in (1) is performed as far as possible, it is easy 
to regain the standard eigenvalues, (3), together with the usual degeneracies. 
For the general operator (2) the ^-function to define is clearly, 

oo ^ 

C(s, a, a; d) = ^ — — — -j , (4) 

and if d = 1, giving the hemisphere, we have. 

From this, it is easy to confirm that the sum of the Dirichlet and Neumann 
functions on the hemisphere does equal the full-sphere ^-function [16]. A technical 
point is that for minimal coupling and Neumann conditions, the origin, m = 0, is 
to be omitted from the sum. 

There are at least three approaches to the continuation of the ^-function (5). 
One is to use a summation formulae, such as Plana's (see [35] for this specific 
case). The second is that used by Minakshisundaram and later by Candelas and 
Weinberg [36] in a treatment of sphere ^-functions and employs a Bessel function 
identity to perform the summation over m. The third method involves an expansion 
in the associated parameter, a^, and, after some manipulation with the binomial 
coefficient, the problem is thrown onto the continuation of a series of Hurwitz 
functions, a standard matter ^. 



^ In this case the dj are the integer degrees associated with the polytope symmetry group, F. For 
the simplest (cyclic) case, Z^, = q with the rest unity. 

^ For a = 0, i.e. the Barnes ^-function, Barnes [37] gives the reduction to a finite sum of Hurwitz 
C-functions and remarks that it could be made fundamental for the theory when d = 1. This 
reduction has been rediscovered in many, more recent discussions of the multiple F- function, 
restricted, as they are, to just this case. In particular, on p. 432, Barnes gives an expression 
for logF„ essentially equivalent to formulae of Vardi [13] and Kanemitsu et al, [38]. A similar 
reduction is also possible for rational degrees. 
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The first approach is not practicable for general degrees, d 7^ 1, (4), since there 
are several genuine summation variables. The second method yields the Bessel 
function expression 

which could be continued to give the analogue of (1) or treated as in Candelas 
and Weinberg [36] and Chodos and Myers [39] to enable values at negative s to be 
computed. 

This is a practical method of obtaining a continuation of the sphere (^-function 
but as our main interest is really with the Barnes function we apply the third 
procedure and expand in o? . One then encounters the continuation of an infinite 
series of Barnes ^-functions, suitable only for special values of s but sufficient for 
our purposes. 



3. Sphere determinants. 

The first step is a standard expansion in powers of o? ^ {cf [40]), 

C(., a, a I d) = ^2rS± + lh^_is + r^ ^^^2, + 2r, a | d), (7) 

r=0 

which allows certain information to be extracted, but which does not constitute a 
complete continuation. The principle being applied here is expansion in terms of 
'known' ^-functions and the continuation of the Barnes ^-function will be assumed 
in principle to be already achieved by Barnes with the numerical computation of 
any particular case looked upon as a soluble technical challenge. 

Barnes has given the values of Cdi^^ ^ I ^) * ^ nonpositive integer, and 
also the residues and remainders at the poles s = 1, . . . , (i, in terms of generalised 
Bernoulli functions rapidly computed by recursion. Hence, for example, one can 
easily find the values C,{—n, a, a | d) for n = 0, —1, . . ., [23]. 

The derivative at also follows, but not so directly. Firstly from (7) 

C'(0, a, a I d) = 2C^(0, a | d) + ^ ^ (r^M + \N2Ad) £ \ 

r=l \ k^l 

00 2r 

+ E ^Cd(2r,a|d). 

r=u+l 
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The notation is that N^{d) is the residue and R^{d) the remainder defined by the 
behaviour at the known Barnes poles^ 

Us + r,a\d)^ ^^^^ + RJd) ass^O, (9) 
s 

where 1 < r < d and u = d/2 if d is even and u = {d — l)/2 if d is odd. 
The problem is the evaluation of the infinite series on the right. 

4. (^—function sums. 

Expression (8) concentrates our attention on the sum of Barnes ^-functions, 

oo 2r 

J2 ^a2r,a|d), (10) 

r 

which can be obtained from 



r 

r=u+l 



E ^ar,a|d), (11) 

r=d+l 



by averaging over ±a. 

For the Hurwitz ^-function {d = 1) this sum is standard {e.g. [41] p. 276) ^ 
but the analysis extends to the more general case as shown by Barnes who has, ([1] 
p.424 eq(l)), (c/[43]), 

^dir, « I d) = log - E TT • (12) 

r=d+l '^^ ' r=l 

The sum on the left is over those ip functions defined by the basic summation 
formula while that on the right contains ip^s that have to be regularised. 

In this paper we maintain Barnes' original notation so that the multiple F- 
and V'-functions are defined by 

C(0,a|d) = logl^, V?^(a) = ^logF», (13) 



^ It may happen that a residue vanishes. 

More complicated sums involving the Hurwitz (^-function have been extensively investigated 
by Srivastava and others, e.g. [42]. Several methods have been adopted and can give different 
'final' expressions leading to identities. 
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where the (d + l)th F-modular form, p^, is given by 

hm [C^(0,a|d)+loga] =-logp,(d). (14) 

Before generahsing (12), we derive it by the method of [23] for a reason ex- 
plained later. 

The integral representation of the Barnes ^-function allows the left-hand sum 
in (12) to be written, using an intermediate regularisation, 

.^..rT{r) Jo ^ ' 

lim r dr [ exp(-aT) - ^ j t'-^K{t) (15) 



r=d+ 



r=0 
d 



= lim (t{s)Us, a + a I d) - ^ ^^r(s + r)Us + r, a | d)\ 



where the 'heat-kernel' K{t) is defined by 

exp(— ar) 



K{t) 



(l - exp(-rf,r)) 

Since the total quantity in (15) is finite, the individual pole terms that arise 
in the s ^ limit must cancel yielding the identity between generalised Bernoulli 
polynomials, 

0(0, a + a I d) - Cd(0, a\d) = Y, t^N^). (16) 



r=l 

The finite remainder is the required result and equals, 

d 



CM a + a\d)- CdiO, a\d)-J2 ^^^r (^) 

- 7(Cd(0, a + a I d) - Cd(0, a \ d)) - J] SZ^^^r)NM 



r=l 



which, in view of (16), yields 
^-^Cdir,a\d) 

r=d+l 



J. ^ J. 

r=l r=l 



d 



(17) 



where is the harmonic number, = Ylk=i 1/^; with Hq = 0. 

To compare with (12) the form of the remainder R^{d) given in Barnes is 

and we see that the two expressions for the sum, (12) and (17), agree. 

We will return to the sphere derivative (8) in section 8, but proceed to generalise 
the sum (12) by differentiating with respect to a, multiplying by and integrating 
back to get, 

E ^ C.(na|d)= / d«,»*(^<"(a + »)-^<"(a)) 

. (18) 

a r+A 

-E(, + A)(r-1)!^^'^(^)- 

This formula generalises one given essentially by Nash and O'Connor [44] for 
the Hurwitz case, d = 1, when it follows after simple geometric summation. (See 
also Choi and Nash [45].) 

Clearly one could continue playing the same game and derive similar summa- 
tions but further progress can be made with (18) when A = n = 0, 1, 2, . . .. We 
concentrate on the integral on the right-hand side which we write in the form 

pa pa 

/ dww'^'ilj''J'\a+w) = a^logT^{a+a)-6^QlogT^{a)-n dw w"^'^ logT ^{a+w) , 
Jo Jo 

(19) 

with the understanding that the final term vanishes when n = 0. 

5. Moments of logF and ip. 

As a preliminary, we discuss integrals of the type 

/•a pa 

/ dww''-UogTa{w), / dt(;t(;"V?^(t(;), (20) 
Jo Jo 

which can be treated by the method in [46] since the algebraic technique applies 

equally well to the Barnes C^function, or indeed to any 'modified' C^function of the 

form,^ 

((g,«;) = y ^ , (21) 



^ This occurs in many places and w could be thought of as a (mass)^ or as a Laplace trans- 
form / resolvent variable . 
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as it relies solely on iteration of the basic relation 

dC{s,w) 



dw 



= -sC(s + !,«;). (22) 



The formulae in [46] give indefinite integrals. We here choose the definite forms, 
(20) with a lower limit of zero. In the Barnes case, the only difference is that log v^27r 
becomes log but this can be avoided formally by extending the summation in [46] 
eq.(109) down to / = and then making some algebraic transformations to arrive 
at the result, 



If" ^ ( —1)^ 

- / dww-i;^^\w) = -J2^—J-{Cd+l-^{-n,d, I d,) + H^C,_,,_,{-n,d, \ dj) 
nl Jo .^^ nl 

(23) 

The first term on the right-hand side of (23) comes from the lower limit, = 0, 
which we have treated by iterating the basic recursion, 

Q{s,w\d) = Cd{s,w + d^ I d) + Cd-i(s,«^| dj, (24) 

where d^ is any degree and is the [d — l)-vector obtained by omitting the d^ 
component from d. Equation (24) has been iterated on d down to Co("^: — 
For convenience we have chosen the d^ to be di,d2,---,d^ taken in turn. The 
notation in (23) then is that the vector d^ has components ((ij, . . . , d^), e.g. d^^ = d 
and is just the single number d^. Finally we let w tend to zero. Incidentally, 
(24) reveals why d^ is a quasi-period.^ 

Adamchik derives the result (23) for the simpler Hurwitz function {d= 1) from 
a series expression for logr(l + x). Espinosa and Moll [47] use recursion to arrive 
at the Hurwitz results. 



6. The general function. 

It is clear from its structure that (12) is the regularised form of an eigenvalue 
sum and its derivation can be paralleled formally for the general (^-function, (21), 

^ For the standard Hurwitz (^-function, the single quasi-period is 1, while for the general case, 
(21), there are none. 
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the precise form of the heat-kernel in (15) not being required. The result is exactly 
equation (17), rewritten in a slightly different notation/-*^ 

f; ^-^ar,w)^log^^^^±^-J^^ (25) 

r=[n] + l ^ ' r=l 

where fi is the order of the infinite set A^, which could be any sequence of numbers, 
e.g. [14], but which most commonly arises as the spectrum of a self-adjoint operator, 
typically the Laplacian on a Riemannian manifold when ^ = d/2, where d is the 
dimension of the manifold. However, one should not limit the meaning of the 
to this case. As they have arisen here, they could be the quantities m.d in (4) and 
might be interpreted as the spectrum of a pseudo-differential operator, such as the 
square root, A-^/^, or something similar. The m.d can be realised as the eigenvalues 
of the pseudo-operator —id. |V| on the d-torus, the degrees, dj, being the inverse 
radii and fx being the dimension, d. In symmetrical cases, for example when several 
degrees coincide, a partially spherical realisation is possible, in accordance with 
section 2. 

The pole residues are given in terms of the coefficients in the short-time ex- 
pansion of the heat-kernel, X] exp ( — {X^ + w)t) (including the mass-squared term, 
w) in a standard way, 

^^w = ^- (26) 

The finite part, FP({r,w), is just another symbol for the remainder, i?^, at the 
possible s = r pole, 

C(s + r,w)^ KM + as s ^ . (27) 

For a given sequence, A^, there may be no relevant poles at s = r G Z so that 
FP({r, w) = C(r, w) and there would be no need to separate the summations as in 
(25). 

The F- functions in (25) are defined in the usual manner, of (14), by 

C'(0,«;) = log3!^, Iogp = -C'(0), (28) 
P 



The identity, (16), becomes a standard one between heat-kernel coefficients. 

An alternative derivation using regularised sums is contained in the next section. 
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with the 'massless' ^-function, denoted by a tilde, 

any zero modes being omitted. With this convention, one could set = C(s,0). 
If '0-functions are also defined in the usual fashion, 

f)<i 

V'^^n«^) = ^iogrH, (29) 

so that for q > fJ,, 

the general formula, (25), looks exactly like the Barnes formula, (12), 



E ^c(^.-)=iog^^-E^^«w, (30) 

leading to the generalisation of (18), for example. 



r=[M]+l 



(31) 



One can also easily produce the formula analogous to (23) for the general 
function, (21), 



- / d^«;"v^a)(^) = -^^(c'(-n) + i?„c(-n)) 



n , (32) 

From the above results one can conclude that the combination C,'{—n,a) + 
H^({—n,a) is significant. The best way of inverting (32) for this quantity is to 
compute the multiple integral 

^—^ / dw{a-wri;('\w)=^—^{C'{-n,a) + H^C{-n,a)) 



n 



-E^^?^(c'(-o+i/,«-o), 

1=0 ^ ' 

(33) 
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where we have adopted the convention^^ that 

dwip^^\w) = logr(Q!) 



so that (23) and (32) hold for n = 0. We confirm this for the Barnes case, equn. 
(23), the right-hand side of which is, at n = 0, 

d 

Cd(0,a|d)-^Cd+i-i(0,rf,|d,). 

i=l 

Using the definitions (13), this can be written 

pMi) r,(j, |di) r,_,id,\d,) T,id,\d^) 



log 



Xdio^ I di) Pd(di) Pd-i(d2) Pi(dd) 
= logrrf(a I di) 

where we have used the fact that T^{d^ \ d) ~ Pd-ii^*) I'lid^ \ d^) = 1. This 
is really a check of algebraic accuracy only, since these relations follow from the 
recursion (24). 



7. Regularised products and sums. 

With (25), or (30), contact has been made with the notion of regularised prod- 
ucts and sums (see e.g. [14,9,10]) because the left-hand side is nothing other than 
the Weierstrass regularisation, 5^* log (l -|- a/{X^ + w)), of the eigenvalue sum 
Elog(l + «/(A^ + «;)), 



(34) 



with 



P{x) = x+ — + ...+ 



2 [^] • 

This amounts to the subtraction of the first [//] terms in the cu-Taylor expansion of 
log (l + a/(A^ + «;)). 



■'^^Adamchik [48] appears to do the same thing in his Proposition 3. 
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The Barnes equation, (12), is a generalised canonical product expression for 
the multiple F-function. 

Equation (8) is, of course, an example of (25) with the eigenvalues of a 
Laplace-type operator. In fact, differentiation of the a-Taylor series of (21), with 
w — > w+a, constitutes another approach to the derivation of (25) (c/ [40] eqn.(9) for 
the Hurwitz function). In the general case we may refer to (25) as Voros' relation, 
[14], eqn (4.12). 

Equation (25) can be obtained by differentiation via a slightly different route, 
[49,43]. First, one defines the Weierstrass a-regularised sum, 



E* 1 
m{X^ + w + ay 



Ell 1 / —s \ a 

\ (\ A_nii-\-r\\s (\ -\-nii\s 



(A^ + «; + «)« {\^ + wy JiK + w)^^' 

(35) 

where sufficient terms in the a-Taylor series have been removed in order to ensure 
convergence for any selected value of s. The integer M determines how many terms 
are to be subtracted. In terms of the order, ^, this means that M = [// — s].^^ 
Differentiation with respect to s, and setting s to zero, produces^'^ 

C\0,w,a)^-y* In (l + a/(A^ + «;)). (36) 

Further, the summation over A^ in (*{s,w + a), (35), can be performed to give 
the continuation 

^ r-s\ 

Cis,w,a) = Cis,w + a)-C{s,w)-J2[ u WC{s + k,w). (37) 

k=l ^ ^ 



In particular 

o H /_ \ 

C*' (0, w, a) = C'(0, w + a)- C'(0, w)-—J2[ , J «'C(s + k, 
which is just (25), after evaluation of the final term. 



w) 



(38) 



s=0 



^"^ This construction is slightly different from that used by Quine et al [9] in an equivalent analysis. 
They set M — [fx] for all s. This is just sufficient to encompass s = 0. Our definition is more 
flexible, as we will see, and is in keeping with the work of Dikii [50] and Watson [51]. 
It is therefore consistent to write C,*' = C,'* . Note also that C,*{s,w,q.) ^ Q:,w). 
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This slicker derivation avoids the use of the heat-kernel form of the ^-function, 
which really acted only as an intermediary. The present approach can be thought 
of as the Mellin transform of the previous one. 

A similar technique, in a particular case, is used by Nash and O'Connor [44] 
App.A. 



8. Sphere determinants again. 



The general sum in (25) was motivated by the example (11), needed in the 
computation of a specific (sphere) determinant. The use made of it depends on 
how much is known about the spectrum. For the direct computation of the sphere 
determinants, obtained by combining (8), (10), (11) and (17), the expression in 
(17) is calculationally explicit since the relevant quantities can be found from the 
properties of the Barnes (^-function. In [23,16] the final expressions were given in 
terms of the derivatives of the Riemann ^-function at negative integers. In the 
sphere combination, the nonlocal Barnes remainders, R^., cancel, in contrast to the 
computation of Casimir energies, for example. 

The further analysis of the hemisphere derivative, (8), can be carried out less 
specifically by starting from 



Z{s,w,a) = J2 ^ % ^ (39) 



m 



instead of from (4). Then 

[m/2] 2r 



Z'(0,«;,a) = 2C'(0,«;)+ ^ ^ (r^,{w)+^-N^,{w)H^_^ 

(40) 



OO 9*, 



+ — C(2r,«;), 

r=[M/2]+l 



where (■(«, w) is the general ^-function, (21). 



The calculation uses the aforementioned reduction of the Barnes function for unit degrees into 

Hurwitz functions. However, what constitutes a 'final' formula is arguable since the Barnes form 
could justifiably be regarded as the end, seeing that what remains is 'merely' a numerical or 
cosmetic affair. For the full-sphere, there is the minor point that one has to compute both the 
Dirichlet and Neumann hemisphere determinants. 
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The final sum in (40) follows from the sum in (25) by averaging, 
E ^C(2r,«;) 

r=[M/2] + l 



and the combination with (40) yields the mentioned cancellation of the ({s, w) 
remainders, and also of the logr^(u'), leaving the formal, but definite, expression 

^//^ ^ , T(w + a)T(w — a) ^sr^ cn'^^ / ^ /.^x 

Z'{0,w,a) = log^ V J2 hO_^NM (42) 

P r=l ^ 

where is the odd harmonic number, = Ylk=o 1/(2^ + !)• 

For the d-hemisphere, this result is given in [23] and, as mentioned there, 
it illustrates the fact that the determinant of a product is not the product of the 
determinants, at least not if the determinant is defined by ^-function regularisation. 
The eigenvalues in (39) factorise, 

(A„ + w)^ - ^ {X^ + w - a) {X^ + w + a) (43) 

and the first term on the right-hand side of (42) gives the product of the determi- 
nants of the individual factors. The remainder is a correction or 'anomaly' which 
was first noticed in physical contexts by Allen [52] and by Chodos and Myers [53] 
and has attracted more recent mathematical, and physical, attention. It is trivially 
zero when a = 0. 

For the hemisphere, = m.d and C(s,'?^) is the Barnes ^-function with unit 
degrees. If the Neumann and Dirichlet hemisphere expressions are added, so as 
to give the full-sphere result, the anomaly contributions cancel in odd dimensions, 
as can be specifically checked. Actually, for odd-dimensional closed manifolds, this 
vanishing is a general result and follows from properties of the heat-kernel expansion 
coefficients. 

Everything is quite definite in (42), which could be taken as the final answer. 
Expressions for the hemisphere Laplacian determinants {a = {d—l)/2) were given in 
[23] in terms of the Barnes function. However, numerical calculation might require 
one to go further and express everything in terms of the Hurwitz or Riemann 
function. (This is arguable.) The details are in [16] where a concrete formula is 



= log 



T{w + a) r{w — a) 



(41) 



r=l 
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produced, 

Csm = 7^ E (1 - + s^'^) a-k) 

'-k=o (44) 

- E ^^""^^^^ + - 1) . 

r—l 

Here, ^^j''^ are Stirling numbers and is the sum of the Dirichlet and Neumann 
hemisphere ^-function residues. 

These residues are given in terms of generahsed Bernoulli polynomials for which 
there exists the handy calculational form, 

d-2 



x={d-l)/2 



and it is easily confirmed that this is zero for odd d, although we actually know this 
at an earlier stage. 

The graph of the results, up to dimension 23, for the Laplacian determinant 
shows a curious difference between odd and even dimensions. The values diverge 
as d increases. This might not be surprising as there are fundamental differences 
between odd and even dimensional spaces, spheres especially. Other values of a in 
(2) can be treated in a like manner. 

We now make some comments on related, full-sphere calculations. Choi and 
Srivastava [18], following Voros, use (25) directly for the two- and three-sphere. 
The Laplacian eigenvalues are written in the usual form, and the standard binomial 
degeneracies are used. Equation (25), with = 0, is then taken as an equation for 
logr(Q;) in terms of logr(O), which is easily found, 

[m] (_ y / \ 

logr(a) = logr(o) + E {Rr + H^-iK) -E* ^ + ^ • ^"^^^ 

The correction terms require and A''^ to be determined from the same a = 
information, which is also relatively straightforward. The most awkward part is 
the evaluation of the last term, i.e. the regularised sum. This is obtained from 
the definition, (34), after using some specific (^-function summations. Although 
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quite workable for small dimensions, this method obscures the general nature of the 
cancellations that must occur. 

In an interesting paper, Quine and Choi [19] express the determinants through 
regularised sums and produce an explicit formula for the (i-sphere using a method 
comparable to our own, [16], and obtain the equivalents of (42) and (44). They do 
not split the mode problem into Dirichlet and Neumann hemispheres. 

Bessel approach to the general (-function. 

An alternative continuation for the ^-function, (39), uses the Bessel identity 
mentioned in section 2. This gives 

Z{s,w,a) = ^j^ dTK{T)l^—j /«_i/2(«t), (46) 

where K{t) is the 'cylinder' heat-kernel 

K{t) = ^e-(^-+^)^. 

m 

The terminology is a reflection of the fact that the squares, (A^ + w)'^, are 
usually the eigenvalues of a second order Laplace-type operator so that K{t) is the 
heat-kernel of a square-root (pseudo)-operator, with different locality properties and 
possible log r terms in its expansion. 

The form (46) has been exploited by Bytsenko and Williams [54] in connection 
with the multiplicative anomaly. 



9. Implicit eigenvalues. 

In the sphere example, the sum, (11), is simply an intermediate calculational 
quantity. It can be made to play a more important role in the evaluation of deter- 
minants especially when the spectrum is not known explicitly. 

The key idea, [49], is to turn (30) around to give a formula, this time for 
logr(w), 

/ \ [m] r 

logr(«;) =^*log 1 + + logr(«; + a)-J2^ i;^^\w) , (47) 

and then to observe, trivially, that the right-hand side has to be independent of 
a. It thus can be calculated at any convenient value. In our work [49,55,43,56] 
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the natural, infinite mass, fimit a — > oo, was selected. The a-dependence has 
to disappear from the asymptotic form and what remains must be logr(t(;). For 
example the last term in (47) can be disregarded since its a-dependence is explicit, 
in finite terms. Note that one is not required to know anything about the nature of 
the V^""^- 

The asymptotic behaviour of \o<gT{w + ct) as a — > cxo has been determined 
([14,57,58], for example) in terms of the heat-kernel coefficients, C^{w\ and contains 
no a-independent terms, apart from the normalisation, log p. Hence we obtain the 
equation 

DetL>-i = - lim FT* (\ + ^ ^ lim ff + T-\ ^^^pY P{-(^/^m) 

(48) 

with the understanding that only the a-independent part of the right-hand side is 
to be retained. The notation has been streamlined a little by dropping explicit 
reference to the parameter w and taking the A^ as the eigenvalues of the generic 
operator D. Keeping only the a-independent part means that the P-polynomial 
bit of (48) is usually irrelevant. 

This equation is particularly valuable when the eigenvalues are given implicitly 
as the roots of some equation, F{a) = 0, since, under certain conditions which are 
often satisfied, the Mittag-Leffier theorem says that F{a) ~ nA(-'^ + ci/A) and the 
asymptotic behaviour of F is often a known episode in special function theory, if 
one is lucky. 

In practice, the method is not necessarily straightforward as the function F 
often occurs after separation of variables and there are remaining summations to 
be dealt with, but it has been applied to d-balls for scalars and higher spin. The 
method also works nicely for bounded Mobius corners [43] , where, one might note, 
the sum (12) occurs with a the parameter in the Robin boundary condition. 

There are several variants of the above technique but all involve the asymptotic 
behaviour of special functions. The spherical cap has been successfully treated by 
Barvinsky et al, [59], using the asymptotics of Legendre functions and the ball 
has been dealt with in various ways that all require Olver-Bessel asymptotics, e.g. 
[60-62]. 
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10. A Kaluza-Klein interpretation. 



The previous development suggests that we make C'{—n^ a) + C{—n, a) the 
subject of equation (33) but rather than the straight derivatives we try to retain a 
Barnesian formulation and rewrite this quantity in terms of a new F-function. 

The best way of doing this is to introduce the ^-function that yields this combi- 
nation as its derivative at 0. This has formal and manipulative benefits. Therefore 
define the ^-function , 

&\s,w) = ^^^^;^C{s-n,w), (49) 

r(s) 

so that 

^(n) (0, yj) = (^'(_n, w) + a-n, w)). (50) 

The new F-function and modulus, T^'^^{w) and p^"\ are defined by 

and 



logp(^) = - lim (c^''y{0,w) + bAogw) 

ui— >0 \ J 

= -C^"^'(0) (52) 

= -^(c'(-^)+^nC(-n)) 

so that T^^^w) = T{w) and p(°) = p. In (52), 6q is the number of zero modes. 
Corresponding ■^-f unctions are defined by 

i^^^'''\w) = -^logT^^\w), (53) 

to complete the formalism. 

Relation (22) gives the recursion on the additional dimensions, 

A^(-)'(0,«;) = -C(--i)'(0,«;). (54) 

The basic idea is to set (^"^^s^w) in place of C(s,w) in formal manipulations 
and then, if desired, return to ({s,w) via (49). 
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Equation (33) now reads in the new notation, 

logr(")(a) = / dw {a - w)''i;^^\w) - V \ \„ logp^ . (55) 

n\ Jo i^-^y- 

An interpretation of (49) is the foUowing. When the are the eigenvalues 
of the Laplacian on a manifold, A^, the derivative, C'(0), is, up to a factor, the 
one-loop effective action of scalar quantum field theory on the 'space-time', A4. If 
a mass is added then one requires ('{0,m?). 

In the same situation, the quantity we have denoted by log (T^'^\w) / p^"^^) 
is, again up to a factor, the efi^ective action density on the Kaluza-Klein manifold, 
M^*^ xM.. (It is a density in the noncompact factor, R^"'.) It therefore can be looked 
upon, when exponentiated, as determining the functional determinant density on 
this manifold. The Kaluza-Klein eigenvalues are + + where k is a 2n- 
dimensional real vector and the integral over k is taken with a certain Att volume 
normalisation to make (49) true. 

The situation with which we are most concerned is when the A^ are the eigen- 
values of a linear pseudo-operator on Ai, as mentioned earlier for the Barnes ^- 
fmiction. The 'Kaluza-Klein' eigenvalues are X^+k.l+w, where k is a real n-vector, 
with < /Cj < oo. The 'Kaluza-Klein' manifold, this time, is MJ^ x Ai. 

In the specific Barnes case, taking Ai to be the (i-torus, the Kaluza-Klein 
manifold x Ai can be thought of as a (d -|- n)-torus with n infinite radii. 



11. Generalised Glaisher-Kinkelin-Bendersky constants. 



The usual Glaisher-Kinkelin-Bendersky constants, Ap, are defined, [8], as the 
AT-independent parts of the sums. 



N 



log A — lim y log m 



m—l 



N 



N — independent 



= lim log I I 



m 



m—l 



AT— independent 



(56) 

which, for p = 0, is Stirling's formula (with the Stirling constant, Aq — ^/2n). 

The original, [63], Glaisher constant, A-^ = exp (l/12 — C^(— 1)), has surfaced 
occasionally in statistical physics, e.g. [64], App.B, as noted by Voros, [14] and the 
Bendersky constants in vacuum quantum field theory, e.g. [65]. 

It is clear from (56) that logA^ is related to the Riemann ^-function, Cr{~p)- 
The exact connection can be obtained by going via the Hurwitz function, noting 
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that, directly from the definition, 



log^^ = hm {C^(-p,N+1)-C'r{-P,1)) 

^ TV— >oo 



(57) 

AT — independent 



and then using, without thought, the standard, large w expansion of C,fi{s, w), [66], 
to arrive at 

log = -Cki-P) - Hp Cr{-p) ■ (58) 

This simple formula is given by Adamchik [48] . 
In our notation, and normalisation, (see (50)), 

iogAp = -p\{-irc'^y{o), 

and, in this way, one can give the usual Glaisher-Kinkelin-Bendersky constants a 
massless Kaluza-Klein determinantal interpretation. 

The F-modular forms, p*-"^-*, (52), are 'new' constants, and g'enerafo erf Glaisher- 
Kinkelin-Bendersky numbers, 

G = o^"^^ 

can be defined in the general situation to mimic the result (58), which holds for the 
Hurwitz ^-function case (when ( is the Riemann ^-function, ). 

In terms of the F-function, the modular form could be defined as the w- 
independent part of logF^") (tu) as — > oo, so that our new constants are. 



logG„ = logp(") = lim logF(")(t(;) 



(59) 

w — independent 



illustrating the formal similarity to (56). 

For the Barnes case, explicit forms follow from (23), 



logGi^)(d) = -^-^ iCd+i-ii-^^ d^ I + H^Cd+l-^i-n, d, \ dj) . 



1=1 
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12. Kaluza-Klein regularised sums. 



The use of the Kahiza-Klein C^^H^, w) in order to facihtate the construction 
of the sphere derivative, (^'{—n, a, a \ d), (see (4)), is contained in [67]. The presen- 
tation of this in [68] is adapted here to a treatment of the general ^-function, (21), 
using regularised sums. 

The combination (49) can be introduced into equation (37) to give 

^ /_ \ 

C^^> {s, w, a) = C^") {s,w + a)- C^^~> {s,w)-J2[k) "'^^"^ + k,w), (60) 
so that 

C(^)*' iO,w,a) = C^^y (0,w + a)-C^^y{0,w)-^^Y. [kj "'^^"^ + ^' ^) • 

'(61) 

Evaluation of the last term, produces immediately the equation corresponding 

to (25), 

± tf c<">(n ») = log - (^(fl<»> + H,_.iv(")) 

r=lfj,]+n+l ^ ' r=l 

^ r(»)(to) ^r\^ * 



r=l 

where the residue and remainder are defined by 

r(")( 
s 

The equality. 



(62) 



(63) 



r=[/i]+n+l 



is shown exactly as in the usual case [n = 0) by summing over and integrating 
over k after the a-expansion of log (l + Q;/(k^ + + w)) , or of log (l + Q:/(k.l + 
+ ^)) depending on the interpretation of the eigenvalues. 
Equation (62) extends (30) and is our final, formal generalisation. 
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If C^") is viewed only as an intermediate quantity one should return to the basic 
(■-function, ({s, w), = C^°Hs, w), and express the Nr^^ and Fli^^ in terms of the A^^. 
and Rj, in order to rewrite the last summation in (62). 

From the definition (49), the poles of (^^'^\s,w) at s = r in the summation 
range, 1 < r < [/i] + n, divide into two sets, those due to r(s — n)/T{s), i.e. r = 
1, . . . , n, and those coming, possibly, from the (-function, i.e. r = n + 1 . . . + n. 

A straightforward calculation produces. 



^i^(4-) + iy,_iAr(-)) 



[n]+n 



E(—aY(r — n — l)\, , , , ^\ 

^ ^ \, ^-{Rr-nH+Hr-n-lNr-nH) , 

r=n+l 

and so equation (62) can be written entirely in terms of the original (-function, 

C{s,w). 

After some shift in the summation variables and multiplication by factors, one 
finds the relation, 

k=[fj.]+l ^ ^ 

= a'"" {C'{-n, w + a) + H^ Ci-n, w + a)) 



n , s 

- E ( ^'(-^' ^) + ^) 



(64) 



fc=0 

{A 



which can of course be proved directly without difficulty. One way is given in [38], 
p. 13, for the simpler Hurwitz (-function {\i= 1). 



13. Conclusion and comments. 



In this work, using the Barnes (-function, we have presented improved and 
generalised expressions related to recent, and not so recent, work concerning sums 
of (-functions which arise incidentally in the computation of sphere determinants. 
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We have concentrated on the full-sphere results but the power of the Barnes 
function shows up when computing quantities on the orbifold factors, S'^/F, where 
r is a polytope symmetry group. We leave these considerations, as well as other 
factorings such as lens spaces, for another time. 

Extensions of our results to forms and spinors, cf [46], is straightforward and 
has topological and possible M-theory applications, [69]. 

A determinant interpretation of the Glaisher-Kinkelin-Bendersky constants 
has been given, based on the Adamchik form, (58). These have been generalised to 
the Barnes (^-function case and also to the more general ^-function, (21). 
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